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Abstract. We estimate by Monte Carlo simulations the configurational entropy of 
TV-steps polygons in the cubic lattice with fixed knot type. By collecting a rich statistics 
of configurations with very large values of N we are able to analyse the asymptotic 
behaviour of the partition function of the problem for different knot types. Our results 
confirm that, in the large N limit, each prime knot is localized in a small region of the 
polygon, regardless of the possible presence of other knots. Each prime knot component 
may slide along the unknotted region contributing to the overall configurational entropy 
with a term proportional to In N. Furthermore, we discover that the mere existence of 
a knot requires a well defined entropic cost that scales exponentially with its minimal 
length. In the case of polygons with composite knots it turns out that the partition 
function can be simply factorized in terms that depend only on prime components with 
an additional combinatorial factor that takes into account the statistical property that 
by interchanging two identical prime knot components in the polygon the corresponding 
set of overall configuration remains unaltered. Finally, the above results allow to 
conjecture a sequence of inequalities for the connective constants of polygons whose 
topology varies within a given family of composite knot types. 
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1. Introduction 

A long, flexible polymer chain in good solvent can be highly self-entangled [TJ [2] and, 
if a ring closure reaction occurs, or if its extremities are hold tied by some device, 
the entanglement can be trapped as a knot [3J Hj. Moreover, because of the excluded 
volume interaction, a knotted molecule cannot change its topological status, without 
breaking and reconnecting back chemical bonds. This for example is the situation one 
encounters in biological systems where special enzymes, called topoisomerases, can pass 
one strand of the double stranded circular DNA through another and knot or unknot the 
molecule, to facilitate elementary cellular processes [HE]. In general, however, there is 
no spontaneous transition between different knotted statuses and in the most common 
experimental situations the topology of the ring does not change in time. Clearly, the 
presence of topological constraints limits the configurational space available to the ring, 
with a consequent reduction of the entropy of the system, compared to the topologically 
unconstrained case [7] . It is then interesting to precisely quantify this entropy loss and 
to determine how it depends on the particular topology (i.e. knot type) considered. 

Unfortunately, most of the theoretical studies performed so far refer to the ensemble 
in which the rings may assume all the topologies. The reason is that polymer rings in 
good solvent can be modelled as self-avoiding polygons (SAPs or simply polygons), 
which are in turn mapped to a magnetic system at its critical point and studied by 
renormalization group techniques [TJ |8j |9] . This approach has led to the well established 
result that the number of Z(N) of iV-steps SAPs grows, for large N, as 

Z{N) ~ Aji N N a ~ 2 (1) 

where the amplitude A and the connective constant \x are non-universal quantities that 
depend on the microscopic features of the chain while a is a universal exponent given 
by a = 2 — dv, where d is the dimensionality of the space and v the metric exponent 
[5]. In d — 3 dimensions, numerical simulations [TU] give for self avoiding loops the 
estimate v ~ 0.587597(7), and consequently a ~ 0.237209(21), in agreement with field 
theoretical results [8] . Since for the subset of SAPs with a given knot type k the above 
mentioned mapping is not valid anymore, there is no field theory argument to establish 
a scaling similar to ([!]) for Z k {N). However, it is reasonable to expect that 

Z k (N) ~ A k ^N a "~ 2 (2) 

where \i k an d ctk are, respectively, the connective constant and the entropic exponent 
of the subset of SAPs with fixed knot type k. For a generic knot type k there is no 
rigorous relation between fi k and \i but in the case of unknotted polygons (i.e. SAPs 
with trivial topology, k = 0) it is possible to prove rigorously that fi$ < n |7J whereas 
numerical estimates of a$ suggests the intriguing identity a ~ a$ [UJ, although results 
presented so far are not sharp enough to rule out completely a possible, although small, 
discrepancy between the two entropic exponents, i.e. a% ~ a. One among the results 
presented here concerns the improvement of the estimate a — ct% and of the ratio A® /A 
(see section |2|), this one performed, to our knowledge, for the first time. 
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Note that Eq. ([I]) and ^ with k = have interesting implications for the 
probability of realizing an unknot in the ensemble with unrestricted topology, P(/){N) = 
Z 9 (N)/Z(N). Indeed, from Q and @ with k = one gets 

and since /i0 < fi we get the well known result that the unknotting probability goes to 
zero exponentially fast with N [7]. The parameter N = 1/ log(/i//i0) gives a typical 
number of steps above which the unknot probability is reasonably low or, in other 
words, the occurrence of knots is not negligible anymore. Previous numerical estimates 
for polygons on the cubic lattice gave jV f»2 x 10 5 [321 [13], [TJ] . 

Since for polymer rings with a generic, fixed knot type k neither analytical tools 
nor rigorous arguments are available, one has to rely entirely on numerical approaches 
and scaling arguments in the analysis of the above issues. By using the BFACF 
algorithm [TSJ [TO] (the acronym comes from the initials of the authors) coupled to a 
multiple Markov chain sampling technique, and assuming for SAPs with fixed knot type 
k the scaling ([2]), evidence found JT7J that 

l^k = H a k = a 9 + Tr k , (4) 

where ir k is the number of prime components in the knot decomposition of k. It is 
interesting to notice that results similar to ^ have been obtained also for off-lattice 
models of rings such as the bead-rod models [18J suggesting that the scaling behaviour 
@ with Q is a universal property of loops in free space with a given knot type 
k. Relations in Q are consistent with recent findings showing that prime knots in 
swollen rings are weakly localized, i.e. have an average "length" (I) ~ A* with an 
exponent < t < 1, which has been estimated in [191 12Q1 [21] as t ~ 0.7. Indeed weak 
localization of prime knots implies that, in the limit N — > oo, each prime component 
behaves essentially as a decorating vertex fluctuating along the unknotted ring. This 
additional configurational degree of freedom brings a factor N in front of Z$ for each 
prime component and, in the general case of a knot k made by n k prime components, 
one may guess: 

Z k (N) ~N^Z $ (N). (5) 

Although the above simple argument furnishes a plausible explanation of relations Q, 
it is too crude to fully characterize the entropy of a knotted ring even in the large N 
limit. For example, the amplitude A k is still undetermined and there is no trace of the 
type of prime knots that contribute to Z k (N). In fact, by regarding prime knots as 
point-like objects, we are neglecting the effective entropic cost that the system has to 
pay in order to tie them into unknotted loops. This entropic cost would decrease the 
Z k (N) in ([5]) by a factor, say, C k giving the more precise expression 

Z k (N) ~ ^^N a «- 2+ *K (6) 
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It is interesting to notice that, if Ck is related to a sort of entropic cost to pay in order 
to tie a given knot, its value should depend on the knot type k and not only on the 
number of its prime components. If this is the case, Eq. ^ would furnish a more 
fundamental description of the large N behavior of the entropy of a knot since it would 
distinguish the type of knot hosted by the polygon. This description should depend also 
on topological invariants of k other than 7tk- 

It is important to stress that a numerical check of the validity of ^ and, in 
particular, a numerical estimate of Ck as a function of k, is a quite hard task to perform 
because it requires a good statistics of polygon configurations with very large values 
of N. This is particularly crucial for SAPs on discrete lattices, for which a reasonable 
amount of knotted configurations can be sampled only for iV > iVo ~ 10 5 . This is 
probably the reason why no attempts have been made so far to look in more details at 
the asymptotic form (|6]). In this paper we explore this issue by sampling polygons on 
the cubic lattice with N up to 200000. Unlike in previous Monte Carlo's, where the 
sampling was performed in the fixed knot ensemble using BFACF [151 EE] algorithm, 
we decided to sample in the free topology ensemble by using the very efficient two- 
pivot-points algorithm [22] and subsequently to partition the sampled configurations 
according to their topology 

In section [2] we describe the algorithm that we use to sample knotted SAPs and 
the procedure designed to detect knots out of configurations that, for large values of N, 
turn out to be highly intricated. As a first outcome of this investigation we will give a 
sharper estimate both of the difference a — «o and of the ratio A*./ A. This will establish 
a more detailed relation between the subclass of unknotted rings and the full class of 
rings with unrestricted topology. In section [3] we test the validity of ^ and estimate 
Ck as a function of k. This is the main result of the paper: it will be first established for 



the simplest case of prime knots and later generalized to composite knots. Section 3.3 



also includes further conjectures on the connective constants of SAP ensembles with 
restricted topology. Section [| is devoted to discussion and conclusions. 



2. Model, Monte Carlo method and knot detection procedure 

To model polymer rings with excluded volume interaction we consider iV-step SAPs on 
the cubic lattice, i.e. self-avoiding walks with the two extremities that are one lattice 
distant apart. These polygons are sampled in free space by using the two-pivot moves, 
a fixed- N algorithm that has been proved to be ergodic in the class of all polygons 
and shown to be very efficient in sampling uncorrelated configurations |22j. With this 
procedure we generate configurations with N up to 200000. As an example, in Fig. [T] 
we plot a configuration with N = 50000, together with a closer view of part of it. 

Since the pivot moves can change the knot type of polygons, the topology of 
each configuration must be detected by means of some topological invariant. This is 
indeed the most problematic part of the whole investigation since, even in good solvent 
conditions, very long polygons may assume a very intricated spatial arrangement. This 




Figure 1. Equilibrium configuration of SAP on the cubic lattice with N 
steps (above) and a detail of its central part (below). 



= 50000 
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Figure 2. By applying the step-reduction algorithm, based on the BFACF algorithm, 
that preserves the topology the configuration of Fig. [T] is simplified to the one 
represented in top panel and having N = 82 steps. A further simplification of the 
Dowker code allows to identify the knot as the composite knot 3i#5i knot whose 
minimal diagram representation is shown on the right. 



"geometrical" entanglement gives rise, in general, to knot projections with a very large 
number of unessential crossings (from the topological point of view) that severely hinders 
the knot detection algorithm based on the calculation of polynomial invariants |23j . 

To circumvent this difficulty, we simplify each sampled configuration before 
performing its planar projection. This is achieved by applying to the polygon a 
smoothing algorithm that progressively reduces the length of the chain while keeping its 
knot type unaltered (for a similar procedure, see [231 ESI ESJ) - This procedure is based on 
the N- varying BFACF algorithm [T51IT6] and has the nice feature of being ergodic within 
each knot type. We set a sufficiently small step fugacity (i.e. the parameter conjugate 
to N), such that the algorithm induces a rapid reduction in the number of steps of the 
polygon. This simplification technique can reduce dramatically the number of crossings 
encountered in an arbitrary projection. An example of how efficient this simplification 
procedure can be is shown in Fig. [T] where a configuration of initially N = 50000 steps 
is shrunk down to the N = 82 steps configuration of Fig. [2] A further reduction is 
achieved by performing 500 projections and choosing the projection with the minimal 
number of crossings. The resulting knot diagram is encoded in terms of the Dowker 
code [27]. A further simplification of the Dowker code based on Reidemeister-like moves 
is performed. Finally, a factorization of the Dowker code is attempted. This procedure, 
whenever successful, splits composite knots into their prime components. From each 
component of the original Dowker code we extract, by using Knotfind [28], the knot 
type of the original configuration (see Fig. [2] for the example given in Fig. [I]). In this 
way we have been able to distinguish composite knots with up to 5 prime components, 
and with each component having crossing number up to 11 [23 EE]- The unbiased 
sampling with unconstrained topology allows us to estimate the probability Pk = Z^/Z 
of occurrence of a given knot type k and to estimate its configurational entropy with 




respect to unknotted polygons, i.e. the ratio Zk/Z$. Since the statistics of unknotted 
polygons will be used extensively as a reference, it is convenient to start by performing a 
good estimate of Z$(N). This can be achieved by looking at the scaling of the unknotting 
probability ^ as N increases. 

In Fig. [3] we plot In P$ as a function of N. The two lines correspond to two different 
fits of the data. To estimate the difference a - «o we first perform a nonlinear fitting 
(dashed line) of the form a — N/N + 6 In N. This yields a — a® = b = —3 x 1CT 5 ~ 
confirming the conjecture a = «o- If we now assume a$ = a we can perform a linear fit 
(solid line) a-N/N Q . This gives JV = 210400±1300 and a = 0.003(2). The estimate of a 
strongly suggests that within error bars A® = A. This last result is quite interesting since 
it strengthens the relation between the statistics of unknotted SAPs and the one of all 
SAPs, not only at the level of the entropic exponents, but also at the level of amplitudes. 
Clearly the main difference relies on the entropies per monomer /1 and However, 
the difference \i — \i§ ~ h/Nq is very small: with the most recent and precise estimate 
fi = 4.684044 ±0.000011 by Slade and coworkers [22], we estimate fj, - // = 0.000022(2) 
(that is twice the statistical error for /x) and thus \i$ = 4.684022 ± 0.000013. 

3. The entropic cost of a knot 

3.1. Prime knots 

To estimate the entropic cost Ck for polygons with fixed knot type k, we compute 
the ratios Z^{N) j 'Zq(N). Indeed, by assuming the scaling form ^ we expect 
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Zk(N) / Zqi(N) ~ N nk /Ck- Let us consider first knotted SAPs where k is a prime knot. 
Fig. [4] shows the iV-behavior (in log- log scale) of the ratio Z^N)/ Z$(N) for prime knots 
up to 6-crossings. As expected from ([6]), no exponential behavior is observed and the 
scaling ~ N/Ck is confirmed (note that i^k = 1 since we are considering prime knots), 
with a Ck whose value increases as the knot complexity increases. The estimates of 
ctfc — a!0, reported in Table [T] in the second column, are all consistent with the relation 
ttfc — «0 = 1. The estimates of Ck are reported in the third column of Table [T] The most 
striking feature to notice is the simple relation observed between the value of Ck and 
the knot type k. Indeed, from column 4, it turns out that, to a good approximation, 
the entropy cost necessary to host a prime knot k goes like 

C k * 4 k/3 (7) 

where 4 (see last column of Table [T]) is the minimal length required to tie a knot k on 
the cubic lattice [30]. Thus, the entropic cost Ck is intimately related to a "microscopic" 
property of the knot k, that is, the length of its "ideal" representation [3TJ [32] on the 
cubic lattice [3D] . 

It is tempting to interpret Vk = 4/3 as an equivalent number of monomers "lost" 

by the polygon in order to form the knot. For example, the partition function of a 

N—V3 

trefoil, Z 3l , would be described by the exponential factor /i 1 with = 8. In other 
words, a iV-step polygon with a 3i knot has the same configurational entropy (in the 
limit of large N) as an unknotted polygon with N — 8 steps endowed with a sliding 
decorating vertex (the knot). These findings suggest that it is sufficient to know the 
length Ik of a given prime knot in its ideal lattice representation in order to make an 
accurate prediction of its frequency along a swollen ring. 

The factor of 1/3 is quite intriguing and we have no explanation for that so far. 
Clearly it will be important to test further this value by looking at more complicated 
knots. This would require a much larger statistics and consequently much larger values 
of N. Another interesting issue would be to see if the relation is model-dependent by 
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1.002(7) 


227800 ± 1400 


7.989(4) 
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1.23(25) 


60(24) x 10 7 


13.1(2) 


40 
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1.22(13) 


38(12) x 10 7 


12.8(2) 


40 


63 


1.08(22) 


61(19) x 10 7 


13.1(2) 


40 



Table 1. Estimates of the difference a}. — ag (second column) and of the entropic cost 
Ck (third column) for the simplest prime knots. The third column suggest a simple 
relation between the entropic cost Ck and the the minimal length £k necessary to tight 
the prime knots of the fist column as estimated in |30| . 




looking for example at polygons embedded on different lattices. 



3.2. Composite knots 

We now extend the analysis of to composite knots, namely knots made by connecting 
prime component knots (such as the 3i#5i in Fig. [2]). In the most general case we may 
assume k to be composed by the prime knots k±, k 2 , ■ ■ ■ , k m , each appearing respectively 
7Tx, 7r 2 and 7r m times. The number 7Tj represents somehow the degree of degeneracy of 
the knot prime ki in the composite knot k. For the composite knot k the entropic cost 
Ck could, in principle, depend on the set {ki} in a quite complicated way. However, if 
we still assume that, in the large N limit, each prime knot localizes along the chain, 
regardless of the presence of other knot components, we can make the working hypothesis 
that the cost of a composite knot k\j^k 2 factorizes as Ck x #k 2 — Cfci x Cfc 2 (f° r &i 7^ fo)- 
This will bring to the conjecture that, in the N — > oo limit, 

Z k (N) ~ Z $ (N) 

The presence of the factorial terms 1/ (tt*)! can be explained as follows: if in a polygon 
with a given knot type there are a number 7Tj of the same prime knot ki, then if we 
permute the position of these prime knots in the polygon the set of corresponding 
overall configurations remains the same. Since in the original statistics we do not take 
care of this over-counting we have to divide the original partition function by 1/(^)1, 
and this must be done for any prime knot component in the knot decomposition. 
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Figure 5. Partition function of knots involving copies of 3i divided by that of the 
unknot and multiplied by the factorial of the number of 3i components, in log-log 
scale as a function of the chain length. Straight lines are power-law fits (exponents are 
shown close to them). 



With these notations, the full cost of a composite knot is then 

Ck = U^ C k7 (9) 

i 

Eqs. ([§])- suggests that, if we knew the entropic cost necessary to tie of each 
prime component fcj, the number of configurations Zk(N) of the composite knot k could 
be easily deduced, in the large N limit, by looking at the partition function Z$(N) of 
unknotted polygons of the same length. 

We first check Eq. ([8j for composite knots including only copies of the trefoil knot, 
for which we have good statistics up to four prime components. From Eq. Q we expect 
the following relation to hold: 

Z 3 JZ $ ^N/C 3l 

Z 3l#3 jZ % x 2 ~ (N/C 3l ) 2 (10) 

Z 3l#3l#3l /Z, x 3! ~ (N/C 3l f 

^3i#3i#3i#3i jZ % x 4! ~ (N/C 3i y 

In Fig. [5] we show these ratios times the suitable factorials, in log-log scale as a function 
of N. The four straight lines are power-law fits whose exponents agree within error bars 



with Eqs. (10). Moreover, as expected, all fits cross each other at a single point (C 3l , 1) 
with abscissa C 3l ~ 227000. Hence the starting assumption that the total entropic cost 
to tie a composite knot of 7r, prime knots simply factorizes (see Eq. [9]) is crisply verified, 



at least for trefoil knots. Note that, by extrapolating the data of Fig [5] to larger values 
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Figure 6. Partition function of knots including a prime component 4i divided by 
that of the same knots without the 4i, in log-log scale vs N. The straight line is a 
guide to the eye, scaling ~ N. 



of TV it is clear that for N > C 3l it is more convenient, entropically, to tie composite 
knots made by trefoils than forming unknotted polygons. 

The statistics collected for multiple copies of the next simplest knots, like the 4i, 
is not sufficient to repeat the analysis of equation Eq. (10). We can however look 



at Eq. (|8j) in the case in which, in addition to multiple 3i's, other prime knots are 
present. In particular in Fig. [6] we plot Z Ax /Z^ Z 3l#A jZ 3l , and Z 3l#3l#4l /Z 3l#3l . As 
expected, all are consistent with the presence of the term N/C Al in the scaling, with 
C Al ~ 5 x 10 6 . This result extends to components of different knot type the hypothesis 
of entropic independence between prime components in the statistics of polygons with 
a composite knot type k. 



3. 3. On the connective constant of a class of composite knots 

So far the only results available on the limiting entropy of knotted polygons are the 
rigorous inequality \i% < fi and the conjectured identity \i% = where refers to the 
connective constant of the subset of polygons having a given knot type k. The essential 
difference between fj, and \x$ or fij. is that in the first case the (infinite) sum over all 
topologies is taken into account while for \i$ and the topology is kept fixed. 

By exploiting Eq. (j§| it is tempting to interpolate between the extreme cases /i 
and fik by looking at the statistics of particular subsets of polygons in which an infinite 
(although partial) sum over topologies is considered. 

Suppose for example to consider the set of all polygons that can have an arbitrary 
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number of trefoil components tight in: 

Z {il #M*0 = Z '^N) + Z 3l#3l (N) + Z 3l#3l#3l (N) + ■■■ (11) 

In the limit of large N, using Eq. ([8]), we get 

00 1 / N \ n 

Z (3l#) ~(iV) ~ Y,Zt(N)— ^ ~ Z $ (N)e N ^. 

By rewriting the exponential factor as (/i0e 1//c ' 3 i) Ar = ( y ue -1 / Afo+1 / c ' 3 i) 7V = 
{H( 3l #)°°) N , since C 3l > N , we get /x( 3l #)°o > /i0. It is interesting to notice that, if 
we apply the same argument to the set of composite knots made only by 4i knots, since 
C& x > C 3l > No, we will get /i( 3l #)<x, > /i^^oo > In general we would expect that 
given two prime knots k' and k" with Cyi > Cy 

> /i(/c"#)°° > A*0- (12) 

This can be explained by arguing that each prime knot, being localized, brings the same 
entropic gain ~ N, but the simplest ones require less entropic cost to be formed. On 
the other hand the statistics of topologically unconstrained polygons are, in the large N 
limit, dominated by extremely complex composite knots made by an arbitrary number 
of different prime components. It is then interesting to look at a more complex subsets of 
polygons whose topology is characterized by an arbitrary number of 3is and 4is. Clearly 
2'(3 1 #)oo,(4 1 #)°o(iV r ) > Z(3 1 #) 00 (A^) + Z( 3l #)oo(A^) and by applying the same argument we 
obtain //^^oc^^oo = /j^e 1 ^ 03 ^ 1 ^ ^ . Hence in general we should expect a sequence of 
the kind 

A*0 = V3i — /^i = • • • < • • • < /^(5i#)°° < A t (4 1 #)°° < 

< V (13) 

4. Conclusions 

By sampling polygons with iV up to 200000 we have been able to get accurate estimates 
of the large iV behaviour of the configurational entropy of SAPs with a fixed knot 
type k. We have corroborated the belief that in good solvent conditions and in the 
large iV limit prime knots are localized within small regions that slide independently 
along the unknotted part of the polygon. The existence of each prime component k 
requires an entropic cost Ck whose dependence on k turns out to be relatively simple 
and intriguingly related to the minimal knot length i^, i.e the minimal number of steps 
necessary to build a knot of type k on the cubic lattice. The above findings allow 
to write down a general formula for the partition function of an arbitrary complex 
composite knots and to conjecture a sequence of inequalities relating the connective 
constants of polygons with different topologies, including families of composite knots. 
In the future it would be nice to explore more broadly the asymptotic relation ^ and in 
particular to test the robustness of the relation Ck oc exp(£fc/3) with respect to different 
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polymer models. In particular it would be interesting to test it in the case of off-lattice 
polymers where should be replaced the length of the knot in its ideal representation 
conformations [311 E2]. Finally we hope that, inspired by the results presented above, 



the set of conjectured inequalities in (13) could be put on rigorous basis by following 
new approaches to the problem. 
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